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Abstract 

A general argument leading from the formula for currents through 
an open noninteracting mesoscopic system given by the theory of non- 
equilibrium steady states (NESS) to the Landauer-Biittiker formula 
is pointed out. Time reversal symmetry is not assumed. As a conse- 
quence it follows that, as far as the system has a nontrivial scattering 
theory and the reservoirs have different temperatures and/or chemical 
potentials, the entropy production is strictly positive. 
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1 Introduction 



Landauer-Biittiker type formulas i.e. expressions relating the (charge, energy 
etc) currents through mesoscopic systems connected with electron reservoirs 
to the corresponding transmission coefficients have been proved to be a key 
tool for analyzing the quantum conductance in nanostructures. Obtained 
initially for the stationary case by phenomenological arguments jH], [B], |15j 
they have been widely extended and used. As for the derivation, one usually 
assumes that the reservoirs have a lead geometry and in order to make use of 
the asymptotic form of the scattered state the current is evaluated far away 
from the scatterer [7j, [TTj, [TBI, a procedure justified (at least in the sta- 
tionary regime) by charge conservation. However this approach may become 
problematic for other reservoir geometries when the leads are short or even 
inexistent (see e.g. ^3], [0]) or in non-stationary regime. 

At a more basic level one starts from a non-equilibrium statistical me- 
chanics formulation (e.g. linear response theory, NESS theory etc) and the 
problem of proving the Landauer-Biittiker formula is to show that the ob- 
tained formula for the current can be cast in a form in which the structure of 
the mesoscopic system enters only via its transition matrix of the associated 
scattering problem as suggested by the phenomenological derivation. In this 
note we shall consider the formula for the current as given by NESS theory 
(see e.g [22], [12], [2], PH] and references therein). One starts at t = from 
an equilibrium state of the decoupled system (i.e. no coupling between the 
mesoscopic system and the reservoirs) with reservoirs having different tem- 
peratures and/or chemical potentials. At the one particle level the system is 
described by 

N 

H = H s + ^2H j -n = n s ©f=i Hj (i.i) 

where Tis, H$ are the Hilbert space and hamiltonian of the mesoscopic 
system and 7ij, Hj are the Hilbert space and the hamiltonian of the jth 
reservoir. The coupling, described at the one particle level by V, is switched 
on suddenly at t — 0. In the limit t — > oo the system settles down to a non- 
equilibrium stationary state. The currents out from the reservoirs are defined 
as minus the time variation of their charge. Since the electrons are considered 
independent the second quantization machinery allows to write the currents 
as given by the general NESS theory in terms of one-particle objects. More 
precisely if (3j, Hj are the temperature and the chemical potential respectively 
of the jth reservoir in the initial state then the current out from the kth 
reservoir in the "final" steady state is [2j (the charge of the electron is — e): 

Jk = ieTr H (n + n F n n* + [v,n k ]) = ieTr H (n F n ni[v,n k ]n + ) (1.2) 
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where IL, are the orthogonal projections onto Hj in H, n = X^=i 11^-, 

Q + = s- lim jt(H Q +v) e -itHo Uo ( 13 ) 

(we follow the notation in the physical literature and j2Ij) and 

TV 

Here f^^.{x) are the usual Fermi-Dirac functions: 

o*) (i • ^* ) '• (i-5) 

Two remarks are in order here. Firstly, we would like to stress the fact that 
(jl.2j) gives only the steady part of the current; a (quasi) periodic component 
might also exists if H = H + V has bound states (see e.g. [IB] Section 5.5 
for a discussion about this point). Secondly, due to the fact that there is no 
interaction between electrons, ()1.2|) can be obtained by means of elementary 
second quantization theory without the use of the heavy artilery of NESS 
theory. 

In the framework of NESS theory and under appropriate technical con- 
ditions the case dim Hs = 1 (single state quantum dot) has been thor- 
oughly studied in [2|, in particular the Landauer-Biittiker formula has 
been proved and strict positivity of the entropy production established. Ac- 
tually in this case the model is nothing but the well known exactly solvable 
Wigner-Weisskopf Atom model (known also under the name of Friedrichs 
model) for which the M0ller operator as well as the scattering matrix can 
be explicitly written down so that the Landauer-Biittiker formula can be 
directly verified. Moreover an extension of the analysis in [2], JH| to more 
general models similar to the one considered in the present note will be given 
(C.-A. Pillet: private communication) in [3] announced in |2j. 

An alternative way (and sometimes more satisfactory from the physical 
point of view) of computing currents in non-equilibrium statistical mechanics 
is to start at t = — oo with reservoirs at the same temperature and chemical 
potential and with an equilibrium state of the coupled system and then switch 
on adiabatically the bias in chemical potential and/or temperature. Unfor- 
tunately due to the fact that in this case the "perturbation" is not localized 
the problem is much more difficult and it has been worked out only at the 
linear response theory level. In this context the Landauer-Biittiker formula 
has been shown to hold true at the heuristic level by Baranger and Stone [7j 
and rigorously proved for a tight-binding model for reservoirs by Cornean, 
Jensen and Moldoveanu [Til] . 
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Coming back to the formula (jl.2|) the problem is that the M0ller operator, 
Q + , involves only "half" of the evolution from — oo to oo encoded in the 
scattering matrix so one has to show that one can rewrite the current only 
in terms of scattering matrix and the initial equilibrium state. In the related 
context of adiabatic quantum pumps theory it has been proved in [Sj that 
this is indeed the case for the lead geometry of the reservoirs. 

The aim of this note (which is a revision and extension of [20]) is a nar- 
row one: to outline a general argument leading from (jl.2j) to the Landauer- 
Biittiker formula. The argument is entirely elementary and very general: 
it works for an arbitrary geometry of the reservoirs (e.g. half spaces, semi 
infinite leads with arbitrary section etc) and arbitrary mesoscopic systems 
of finite size. Also we allow the reservoirs to be coupled both via the meso- 
scopic system and by direct contacts [12J. Actually, the only thing which is 
needed is a good stationary scattering theory for the pair (Ho, H + V). This 
is consistent with the generality of the phenomenological arguments leading 
to Landauer-Biittiker formula. For the model at hand the strict positivity 
of the entropy (which is a central issue in NESS theory and has been estab- 
lished under various conditions for by far more general models see e.g. 0, 
, P2] and references therein) follows from the Landauer-Biittiker formula; 
while for systems with time reversal symmetry this is straightforward by the 
argument in [2] in the general case the argument is a bit more involved. 

The content of the note is as follows. In Section 2 we specify the model 
and give the needed formulas from stationary scattering theory. Section 
3 contains the argument leading from the formula ()1.2)) to the Landauer- 
Biittiker formula. In order not to burden the simplicity of the argument in 
technical and notational details we shall give it at the formal level and in 
a simple context: two reservoirs with the same simple (single channel) ab- 
solutely continuous spectrum, a C [0,oo), and mesoscopic systems with a 
finite number of states (i.e. dim Hs = M < oo). For a mathematical sub- 
stantiation for various concrete models one has to make precise the technical 
conditions on Hq and V and then check that one can apply the results of the 
rigorous stationary scattering theory as developed e.g. in [I], [21], [23 (for 
stationary scattering theory at the formal level we send the reader to [TB*]). 
In Section 4 we give some straightforward extensions as well as the argument 
for strict positivity of the entropy production. 
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2 The model and its scattering theory 

As already said in the introduction, at the one particle level the uncoupled 
system is described by 

N 

H = H S + Y, Hj- H = H S ©f =1 Hj (2.1) 

3=1 

where Hs, H$ are the Hilbert space and hamiltonian of the mesoscopic 
system and Hj, Hj are the Hilbert space and the hamiltonian of the jth 
reservoir. We suppose that the spectral representation of Hj j = 1,2 is 
given in terms of generalized eigenfunctions, \ip°j e >> living in an appropriate 
"weighted" Hilbert space, /C* (Gelfand triplets structure: JCj C Hj C /C*) : 

Hj\^l E >= E\<iIjI e >, E G a(Hj) =a C [0,oo). (2.2) 

For / G Hj we denote by f{E) its generalized Fourier transform: 

f(E)=<^ j>E J>. (2.3) 

As concerning V, we suppose to have the following structure (in the 
decomposition given by (jl.ljO : 

/ Vsi V S2 \ 

V = V£ x V 12 . (2.4) 

KVn v* 2 ) 

Since Hs has finite dimension, V$\ and V$2 are finite rank operators. We 
suppose the "direct contact", Vi 2 , to be also of finite rank. Accordingly: 

rrij<M 

Vsj= ^2 v oi\ s ii >< fill ( 2 - 5 ) 

m<oo 

Vi 2 = ^2 v i\9u >< 92i\ (2.6) 
i=i 

where j — 1,2 , {fji}^, {Oj^JLv { s ji}i=i are orthonormal systems in Hj 
and Hs respectively and Vji, vi > 0. 

Since V is of finite rank, by Kato-Kuroda-Birman theory [21], [21], [2H] 
the M0ller operators 

Q± = 8- lim e it{Ho+v) e- itHo U (2.7) 
5 



exist and are unitary from n 7i onto the absolutely continuous subspace, 

u ac n,oiH. 

We impose further conditions on fji , gji in order to assure that Q± provide 
spectral representations for H restricted to H a cH ie for all E G a ac (H) = 
\Jj =1 cr(Hj), with a possible exception of a discrete set, £, fl± have bounded 
extensions in the orthogonal sum of K,* and 

fo^>=n±|i# B > (2.8) 

are generalized eigenfunctions for H: 

H\^ E >= E\^ E > . (2.9) 

A sufficient condition (which at the price of more technicalities can be weak- 
ened) in the case when Hj are discrete or continuous Laplaceans supple- 
mented with boundary conditions is that fji, gji are exponentially localized 
in space. This condition also implies that the generalized Fourier coefficients, 
fji(E), gji(E) (see 1)2.3)1 ) of fji, g^ are smooth functions of E, a fact which is 
needed in order to apply the principal value formula during the proof below. 
The generalized eigenfunctions satisfy the Lippmann-Schwinger [23] . 
[23] equation: 

>= Ke >-(Ho~Et iOr'Vl^z > • (2.10) 
Consider now the scattering operator 

S = tt*_Q+ (2.11) 
and the corresponding transition operator, T, defined by 

S=l-2i%T. (2.12) 
Since S (and then T) commutes with Hq it has a spectral representation: 

S= ! S(E)dE, T= J T(E)dE (2.13) 

where S(E) is a unitary two by two matrix (we are considering the case of 
two reservoirs with simple spectrum). From the unitarity of S(E) it follows 
that T(E) satisfies the so called optical theorem: 

T(E) - T{E) = -2iriT(E)T*(E). (2.14) 

The basic result of the stationary scattering theory is the formula for T(E) 
in terms of the generalized eigenfunctions of H [TH] , [I] , [21] , 



T jk (E) =< Vifr+s >=< tf^, Vn + i>l E > . (2.15) 
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3 Landauer-Biittiker formula 

To prove the Landauer-Biittiker formula in the context described above 
amounts to prove : 

Proposition 1. 

h = 2eTr w (fi + n F n fi;[y,n 1 ]) = 

-2e^ J ao dE{fl^{E) - fZ 2 {E))\T l2 {E)\\ (3.1) 

The second equality in ()3.1|) is the main result of this note. 
We compute j\ from (jl.2j) . Inserting the formula for V (see (12. 4 J) . ()2.5|) . 
((222)) an d computing the trace in appropriate bases one gets: 

Y!{L7 vi < gu, n + u F u n* +92l >). (3.2) 

Using the spectral representation of Uo-Follo (see (|1.4j) ) in ()3.2j) or, alterna- 
tively, evaluating directly the trace in the r.h.s. of (|1.2|) in the generalized 
basis of H one gets: 

31 = -2e j ao dE{fl D ^{E)^ < Vn+^ILin+rlfa > + 

f^ 2 (E)Z < VQ + r 2 , E , TLiO+r^ >}. (3.3) 

Let us compute first the coefficient of fp^i.E) in ()3.3|) . Using the Lippmann- 
Schwinger equation (see (j2.10j0 for Q + ip® E and the spectral representation 
of H one has: 

< vn + ^ E ^i E >-< vn+^Ui^^vn+ri^E >= 

Tm-JdE' ^^^ . (3.4) 

Now the important fact is that we need only the imaginary part of (j3.4jl . 
Then using the principal value formula 

-J— = i7r 8(0) + PV-, (3.5) 

to evaluate the integral in ()3.4|) ( since | < -?/>i E ,, VVL+ipi E > | 2 depends 
upon E' only via fu(E') and gu(E') which are smooth by assumption, this 
is legitimate) one obtains that the coefficient of {E) in ()3.3|) is 

S(7\^E)-27r|T n (E)| 2 ). 

Now the use of the optical theorem (j2.14|) leads to the conclusion that the 
coefficient of fp^^E) in (|3.3|) is 7r|Ti 2 (£ ; )| 2 . A similar computation for the 
coefficient of fp 2 %(E) in (jSH (in this case the term linear in T(E) vanishes) 
leads to — 7r|Ti 2 (i?)| 2 and the proof of (j3.1|) is finished. 
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4 Generalizations and strict positivity of the 
entropy production 

We give first some straightforward extensions of the result in previous section. 

i. A similar proof applied to the energy current (see e.g. 0) gives: 

$! = -^(^noFono^^n^on!]) = 

2nJ aQ dE(f^(E) - fZ,( E ))E\T l2 {E)\ 2 - (4-1) 

ii. The condition that cr(Hi) = a(H2) (as sets) is not necessary ; in the 
general case only the energies in the intersection of cr(Hi) with <j(H 2 ) can 
have nontrivial scattering and then contribute to the current. 

iii. The straightforward generalization of (j3.1|) . (|4.1|) to the case of N 
reservoirs is given by: 

j k = ieTr H (Q + U F U Q* + [V,U k ]) = 
-2enj ao dEj:U(fZS E ) ~ fp^( E ))\T kj (E)\*, ( 4 - 2 ) 

$ fe = -tTr n {ft + U F U Q* + [V,U k H U k ]) = 
2nj ao dEJ2f =1 (fZ h ( E ) - f^(E))E\T k3 (E)\\ (4.3) 

Notice that for N > 2 and systems without time reversal symmetry one can 
have \T kj (E)\ 2 ^ \T jk (E)\ 2 . Still EtiJk = 0, £f=i $ fe = as required 
by charge and energy conservation, due to the fact that the unitarity of S 
implies that T(E) is a normal matrix: 

T(E)T*(E) = T*(E)T(E). (4.4) 

iv. If at some energy E the spectra of Hj do not have multiplicity one 
then Tjk(E) become operators Tj k (E) : Tlj(E) — ► H k (E) and \Tj k (E)\ 2 in 
(jOD are to be replaced by Tr H][E) T ]k {E)T* k {E). 

Consider now the entropy production. The entropy production rate (as 
given by heuristic thermodynamic arguments) has been identified as the 
(thermodynamic and t — > oo) limit of the relative entropy of the evolved 
state with respect to the initial state, |2Sj, [T7] . and then its posi- 

tivity follows from Klein's inequality. However, although for finite reservoirs 
and times its strict positivity can be also easily obtained from Klein's in- 
equality it is not clear whether it survives when taking thermodynamic and 
infinite time limits. Accordingly the strict positivity of the entropy produc- 
tion has to be established for each concrete realization of NESS theory and 
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many results are known (see e.g. jl], |2] 5 C3 an d references therein). For the 
model considered in this note the entropy production, cr, writes as (see e.g. 
formula (2.61) in ^2] or formula (6.46) in [2J: 



N 



o = - Pk(®k - (4-5) 



k=l 

where 

N 



Ik = ~~jk = 2tt I dEj2(f^ k ( E ) - f^(E))\T kj (E)\ 2 (4.6) 
is the particle current. From (|4.8jl . (J4.2)) and ()4.6)1 and using the notations: 

/( X ) = (l + 2 ;)- 1 ; X E j = fy(E - fJLj) (4.7) 

one obtains: 

N 



a = -2tt f dEj2 (f(xE,k) ~ f(x E , 3 ))x E , k \T kJ (E)\ 2 . (4i 

J (JO A 1 1 

Using 



j,k=i 



N N 

J2\T l>m (E)\ 2 = J2\ T m,im 2 , l = l,2,...,N (4.9) 

m=l m=l 

which follows from (|4.4j) . one can rewrite ()4.8|) as 

iV 

a = 2tt / d£ ^ f{x E j)(x E> k ~ XE^Tk^E)?- ( 4 -!0) 

We shall say that the model has a nontrivial scattering in the channel (jo, &o) 
if Tj 0iko (E) ^ (on a subset of <r of positive Lebesgue measure). As expected, 
the non-triviality of the scattering implies strict positivity of the entropy 
production: 

Corollary 2. Suppose that the scattering is nontrivial in the channel (jo, ko) 
and in addition (/?,-„, /x,-„) ^ {Pk ,^k ) (i.e. at least (3 jo ^ (3 ko or fj, jo ^ ^ k(l ). 
Then 

a>0. (4.11) 
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Proof. For the time reversal symmetric case i.e. \T^ m (E)\ = \T m j(E)\ 
(|4.11|) follows from a simple argument in [2 : by symmetrizing the sum in 
either ()4.8j) or 1)4.10)1 one obtains: 

r - 

a = n dEj2 (fM - f(x E>k ))(xE,k ~ x Ejj )\T k>j (E)\ 2 (4.12) 
u ° j,k=l 

and then ()4.11j) follows from the fact that f(x) is strictly decreasing. In the 
general case the argument is a bit more involved and, as expected, mimic the 
proof of Klein's inequality. Let F(x) be a primitive of f(x). Since f(x) is 
strictly decreasing, F(x) is strictly concave i.e 

F(x)<F(y)+f(y)(x-y) (4.13) 

and the equality is reached only for x = y. From ()4.13|) and ()4.9)) it follows 
that 

Ef^f^EjKxE^-XE^in^E)^ > £^ =1 (F(x £)fc ) - F(x E ,)\T k ,(E)\* 
= Ef, k=1 F(xE, k K\T k ,(E)\ 2 - \T hk {E)\ 2 ) = 0. (4.14) 

and the equality is attained only if for all pairs (J, k) and all (a.e) E, (xE, k — 
XEj)\Tk,j(E) \ = 0. Then the observation that x E ,k ~ x ej = holds true for 
at most one energy finishes the proof . 
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